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Abstract: We introduce a parametrized notion of genericity for Delaunay triangulations which, 
in particular, implies that the Delaunay simplices of ^-generic point sets are thick. Equipped with 
this notion, we study the stability of Delaunay triangulations under perturbations of the metric 
and of the vertex positions. We quantify the magnitude of the perturbations under which the 
Delaunay triangulation remains unchanged. 
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Stabilite des triangulations de Delaunay 

Resume : On introduit une notion paranietree de genericite pour les triangulations de Delaunay 
qui iniplique en particulier que les siniplexes de Delaunay d'ensembles delta-generiques ne peuvent 
pas etre tres plats. En s'appuyant sur cette notion, on etudie la stabilite des triangulations de 
Delaunay en presence de perturbations de la nietrique et de la position des points. On quantifie 
I'amplitude des perturbations pour lesquelles la triangulation de Delaunay reste invariante. 
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1 Introduction 

One of the central properties of Delaunay complexes, which was demonstrated when they were 
introduced [Del34', is that under a very mild assumption they are embedded, i.e., they define 
a triangulation of Euclidean space. The required assumption is that there are not too many 
cospherical points; the points are "generic". The assumption is not considered limiting because, 
as Delaunay showed, an arbitrarily small afhne perturbation can transform any given point set 
into one that is generic. 

Given the assumption of a generic point set, we are assured that the Delaunay complex 
defines a triangulation, but a couple of issues arise when working with these triangulations. One 
is that the Delaunay triangulation can be highly sensitive to the exact location of the points. For 
example, the Delaunay triangulation of a point set might be different if a coordinate transform 
is first performed using floating point arithmetic. 

Another problem concerns the geometric quality of the simplices in the triangulation. We 
define the thickness of a simplex as a number proportional to the ratio of the smallest altitude to 
the longest edge length of the simplex, and we demonstrate why this is a useful measure of the 
geometric quality of the simplex. For points in the plane, if there is an upper bound on the ratio 
of the radius of a Delaunay ball to the length of the shortest edge of the corresponding triangle, 
then there is a lower bound on the thickness of any Delaunay triangle. However, when there are 
three or more spatial dimensions, the thickness of Delaunay simplices may become arbitrarily 
small in spite of any bound on the circumradius to shortest edge length. 

Both of these issues are shown to be related to points being close to a degenerate (non-generic) 
configuration. We parameterize Delaunay's original definition of genericity, saying that a point 
set P C M™ is (5-generic if every m-simplex in the Delaunay complex has a Delaunay ball that is 
at a distance greater than S to the remaining points in P. We show that a bound on S leads to a 
bound on the thickness of the Delaunay simplices, and also that the Delaunay complex itself is 
stable with respect to perturbations of the points or of the metric, provided the perturbation is 
small enough with respect to 5 in a way that we quantify. In a companion paper |BDG13| . we 
develop a perturbation algorithm to produce (5-generic point sets. 

The stability of Delaunay triangulations has not previously been studied in this way. Related 
work can be found in the context of kinetic data structures [AGG"'"10] or in the context of robust 
computation [BS04) . and in particular, the concept of protection we introduce in Section p] is 
embodied in the guarded insphere predicate which has been employed in a controlled perturbation 
algorithm for 2D Delaunay triangulation |FKMS05] . 

Our interest in the problem of near-degeneracy in Delaunay complexes stems from work on 
triangulating Riemannian manifolds. An established technique is to compute the triangulation 
locally at each point in an approximating Euclidean metric, and then perform manipulations 
to ensure that the local triangulations fit together consistently [BWYlTl IBGllj . The reason 
the manipulations are necessary is exactly the problem of the instability of the Delaunay trian- 
gulation, and sometimes this is most conveniently described as an instability with respect to a 
perturbation of the local Euclidean metric. 

Although we make no explicit reference to Voronoi diagrams, the Delaunay complexes we 
study can be equivalently defined as the nerve of the Voronoi diagram associated with the metric 
under consideration. We provide criteria for ensuring that the Delaunay complex is a triangula- 
tion without explicit requirements on the properties of the Voronoi diagram |ES97| , in contrast 
to a common practice in related work |LLOO| lLS03l ICDR05| IDZM08| [CGT2] . 

After presenting background material in Section [2] we introduce the concept of (5-generic point 
sets for Euclidean Delaunay triangulations in Section [3] We show that Delaunay simplices of 
(5-generic point sets are thick; they satisfy a quality bound. Then in Section|4]we quantify how 6- 
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genericity leads to robustness in the Delaunay triangulation when either the points or the metric 
are perturbed. The primary challenge is bounding the displacement of simplex circumcentres. 
We conclude with some remarks on the construction and application of (5-generic point sets. 

2 Background 

Within the context of the standard ?Ti-dimensional Euclidean space M™, when distances are 

determined by the standard norm, ||-||, we use the following conventions. The distance between a 

point p and a set X C M™, is the infimum of the distances between p and the points of X, and is 

denoted dRm[p,X). We refer to the distance between two points a and b as \\b — a\\ or dRm(a, 6) 

as convenient. A ball BT^m{c;r) — {x\ ||x — c|| < r} is open, and i?Rm(c;r) is its topological 

closure. We will consider other metrics besides the Euclidean one. A generic metric is denoted 

d, and the associated open and closed balls are B{c;r), and B{c;r). Generally, we denote the 

topological closure of a set X by X, the interior by intX, and the boundary by dX. The convex 

hull is denoted conv(X), and the affine hull is aS(X). 

If U and V are vector subspaces of M™, with dimU < dimV, the angle between them is 

defined by 

sin Z{U,V) ~ sup ||u — TTyulJ, (1) 

ueu 

||tt||=i 

where Try is the orthogonal projection onto V. This is the largest principal angle between U and 
V . The angle between affine subspaces K and H is defined as the angle between the corresponding 
parallel vector subspaces. 

2.1 Sampling parameters and perturbations 

The structures of interest will be built from a finite set P C M™, which we consider to be a set of 
sample points. If D d M™ is a bounded set, then P is an e-sample set for D if rfRm(a;, P) < e for 
all X € D. We say that e is a sampling radius for D satisfied by P. If no domain D is specified, 
we say P is an e-sample set if dTRm{x, P U 9conv(P)) < e for all x G conv(P). Equivalently, P is 
an e-sample set if it satisfies a sampling radius e for 

D,{P) = {xe conv(P) |dR™(x,9conv(P)) > e}. 

The set P is X-separated if dRni(p, q) > X for all p,q (z P. We usually assume that the sparsity of 
a e-sample set is proportional to e, thus: A = iioe. 

We consider a perturbation of the points P C M™ given by a function ^ : P — > M™. If C 
is such that (iR™(p, C(p)) ^ Pj '^^ say that C is a p-perturbation. As a notational convenience, 
we frequently define P = C(P)> ^-^id let p represent C,{p) E P. We will only be considering 
p-perturbations where p is less than half the sparsity of P, so ^ : P ^ P is a bijection. 

Points in P which are not on the boundary of conv(P) are interior points of P. 

2.2 Simplices 

Given a set of j + 1 points {pq, . . . ,pj} C P C M™, a (geometric) j -simplex a = [po, • • • ,Pj"] is 
defined by the convex hull: a = conv({po7 • • • jPj})- The points pi are the vertices of a. Any 
subset {pi„, ■ . ■ ,Pii^} of {pq, . . . ,pj} defines a /c-simplex r which we call a face of a. We write 
r < cr if T is a face of a, and r < a if r is a proper face of a, i.e., if the vertices of r are a proper 
subset of the vertices of a. 
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The boundary of a, is the union of its proper faces: da = lJr<(T '''■ ^^ general this is distinct 
from the topological boundary defined above, but we denote it with the same symbol. The 
interior of a is int a = a \ da. Again this is generally different from the topological interior. In 
particular, a 0-simplex p is equal to its interior: it has no boundary. Other geometric properties 
of a include its diameter (its longest edge). A(ct), and its shortest edge, L{a). 

For any vertex p G cr, the face oppposite p is the face determined by the other vertices of a, 
and is denoted ap. If r is a j-simplex, and p is not a vertex of r, we may construct a (j + 1)- 
simplex a — p * t, called the join of p and r. It is the simplex defined by p and the vertices of 
r, i.e., T = ap. 

Our definition of a simplex has made an important departure from standard convention: 
we do not demand that the vertices of a simplex be affinely independent. A j-simplex ct is a 
degenerate simplex if dimaff(CT) < j. If we wish to emphasise that a simplex is a ^'-simplex, we 
write j as a superscript: cr-'; but this always refers to the combinatorial dimension of the simplex. 

A circumscribing ball for a simplex a is any m-dimensional ball that contains the vertices of a 
on its boundary. If a admits a circumscribing ball, then it has a circumcentre , C{a), which is the 
centre of the smallest circumscribing ball for a. The radius of this ball is the circumradius of a, 
denoted R{a). In general a degenerate simplex may not have a circumcentre and circumradius, 
but in the context of the Euclidean Delaunay complexes we will work with, the degenerate 
simplices we may encounter do have these properties. We will make use of the affine space N{a) 
composed of the centres of the balls that circumscribe a. This space is orthogonal to aff (tr) and 
intersects it at the circumcentre of a. Its dimension is to — dim aff (cr). 

The altitude of a vertex p in cr is D{p, a) = (iK"i(p, aff (cTp)). A poorly-shaped simplex can be 
characterized by the existence of a relatively small altitude. The thickness of a j-simplex a is 
the dimensionless quantity 



'^i'^) ~ \ .„:„ D(p.<y) 



ifj^O 



minpgo- j^(ct) otherwise. 

We say that a is To-thick, if T(ct) > Tq. If a is To-thick, then so are all of its faces. Indeed if 
T < a, then the smallest altitude in r cannot be smaller than that of ct, and also A(r) < A(ct). 

Our definition of thickness is essentially the same as that employed by Munkres |Mun68| . 
Munkres defined the thickness of ct-' as ^raT; where r{a^) is the radius of the largest contained 
ball centred at the barycentre. This definition of thickness turns out to be equal to ^irTT(CT-'). 

Whitney |Whi57| employed a volume-based measure of simplex quality, and variations on this, 
typically referred to as fatness, have been popular in works on higher dimensional Delaunay-based 
meshing [CDE +OOllLiOBIlBWYlll . We find a direct bound on the altitudes to be more convenient, 
because it yields a cleaner and tighter connection between the geometry and the linear algebra of 
simplices. Typically, a bound on some geometric displacement related to a simplex is obtained 
by bounding the inverse of a matrix associated with the simplex, and thickness is well suited for 
this task. 

As a motivating example, consider the problem of bounding the angle between the affine hull 
of a simplex and an affine space that lies close to all the vertices of the simplex. Such a bound 
is relevant when meshing submanifolds of Euclidean space, for example, where it is desired that 
the affine hulls of the simplices are in agreement with the nearby tangent spaces of the manifold. 

Whitney jWhi57l p. 127] obtained such a bound, which manifestly depends on the quality 
of the simplex. Using thickness as a quality measure we obtain a sharper result with a simpler 
proof: 
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Lemma 2.1 (Whitney angle bound) Suppose cr is a j-simplex whose vertices all lie within 
a distance 77 from a fc-diniensional affine space, H C M™, with k > j. Then 



smZ{aS{a),H)< ^^ 



T(a)A(a)- 



The idea of the proof is to express the unit vector u in Equation (fl]) in terms of a basis for aff (cr) 
given by the edges of a that emenate from some arbitrarily chosen vertex. The projection u of 
u into H can then be expressed in terms of the projected basis vectors, using the same vector 
of coefficients. Since the vertices of a all lie close to H, the projected basis vectors do not differ 
significantly from the originals, so bounding the magnitude of the difference between u and u 
comes down to bounding the magnitude of the vector of coefficients of the unit vector u. This 
bound depends on how well-conditioned the basis is, and this is closely related to the thickness 
of a. 

These observations can be conveniently expressed and made concrete in terms of the singular 
values of a matrix. An excellent introduction to singular values can be found in the boolc by 
Trefethen and Bau |TB97I Ch. 4 & 5], but for our purposes we are primarily concerned with the 
largest and the smallest singular values, which we now describe. 

We denote the i*'' singular value of a matrix A by Si{A). The singular values are non-negative 
and ordered by magnitude. The largest singular value can be defined as si{A) — supii^ii j^ ll^2;||; 
it is the magnitude of the largest vector in the range of the unit sphere. The fist singular value 
also defines the operator norm: \\A\\ = Si{A). The standard observation that a bound on the 
norms of the columns of A yields a bound on \\A\\ is obtained by a short calculation. 

Lemma 2.2 If 77 > is an upper bound on the norms of the columns of a, m x j matrix A, then 
\\A\\ < VjV- 

We will also be interested in obtaining a lower bound on the smallest singular value which, for 
an 771 X J matrix A with j < m, may be defined as Sj(A) = infiun^]^ ll^2;||. 

From the given definitions, one can verify that if A is an invertible m x m matrix, then 
si{A~^) = Sm{A)~^, but it is convenient to also accommodate non-square matrices, correspond- 
ing to simplices that are not full dimensional. If ^ is a 777, x j matrix of rank j < m, then the 
pseudo-inverse A^ — {A'^ A)^^ A'^ is the unique left inverse of A whose kernel is the orthogonal 
complement of the column space of A. We have the following general observation: 

Lemma 2.3 If A is a 771 x j matrix of rank j < m, then Si{A^) = Sj_i^i{A)^^. 

The columns of A form a basis for the column space of A. The pseudo-inverse can also be 
described in terms of the dual basis. If we denote the columns of A by {ai}, then the i"^ dual 
vector, Wi, is the unique vector in the column space of A such that wjui — 1 and wjuj = if 
i ^ j. Then A' is the j x m matrix whose i row is wj . 

By exploiting a close connection between the altitudes of a simplex and the vectors dual to 
a basis defined by the simplex, we obtain the following key lemma that relates the thickness of 
a simplex to the smallest singular value of an associated matrix: 

Lemma 2.4 (Thickness and singular value) Let cr = [po, . . . ,pj] be a non-degenerate j- 
simplex in M™, with j > 0, and let P be the m x j matrix whose i"^ column is pi ~ po- Then 

s,iP) > V7T(r7)A(a). 

Proof Let wJ be the 7*'^ row of P^ . By the definition of P^ , it follows that Wi belongs to the 
column space of P, and it is orthogonal to all (pii — po) for i' ^ i. Let Ui = Wi/ \\wi\\. By the 
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Figure 1: Choosing po as the origin, the edges emenating from po in a = [po, . . . ,pj] form a basis 
for aff(CT). The proof of Lemma 2.4 demonstrates that the dual basis {wi} consists of vectors 
that are orthogonal to the facets, and with magnitude equal to the inverse of the corresponding 
altitude. 



definition of Wi, we have wj{pi — po) = 1 = \\wi\\ uj {pi — po). By the definition of the altitude 
of a vertex, we have uj{pi — po) — D{pi,(T). Thus \\wi\\ — D{pi,a)^^. Since Si{A^) — Si{A) for 
any matrix A, we have 



i<i<j 



by Lemma [2. 2| Thus Lemma [2 . 3 1 yields 



,(P) > ^ mm D{p,.,a) = ^Tia)Aia). 



D 



The proof of Lemma |2 .4| shows that the pseudoinverse of P has a natural geometric interpre- 
tation in terms of the altit udes of a, and thus the altitudes provide a convenient lower bound 



on Sj{P). By Lemma 2.2 si{P) < \/?A(cr), and thus T(cr) < 



si(P) 



In other words, T(ct)" 



provides a convenient upper bound on the condition number of P. Roughly speaking, thickness 
imparts a kind of stability on the geometric properties of a simplex. This is exactly what is 
required when we want to show that a small change in a simplex will not yield a large change in 
some geometric quantity of interest. For example, we will use Lemma |2.4| in the demonstration 
of Lemma |4.1[ which is the technical lemma related to the stability of the space of circumcentres 
of a simplex. Lemma [2. 4| also facilitates a concise demonstration of Whitney's angle bound: 



Proof of Lemma 
be the vector su 

and let tt : M" - 



2.1 



Suppose (T = [po, ■ ■ ■ ,Pj]- Choose po as the origin of M™, and let U C M'" 

3space defined by aS{a). Let W be the A:-dimensional subspace parallel to H, 

W be the orthogonal projection onto W. 

For any unit vector u E U, smZ{a.S{a), H) — smZ{U,W) < ||u — ttmIJ. Since the vectors 

Vi — iPi—po)i i G {1) ■ • ■ ) i} form a basis for U , we may write u — Pa, where P is the m x j matrix 

whose z*'' column is Vi, and a £ W is the vector of coefficients. Then, defining X ^ P — ttP, we 

get 

||M-7ru|| = ||Xa|| < ||X|| ||a||. 
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VF is at a distance less than ry from H , because po G W and dR™ {pi,H) < 77 for all < i < j. 



It follows that \[vi — TTViW < 2r], and Lemma 2.2 yields 

ll^ll<2V7^. 

Observing that 1 — \\u\\ — \\Pa\\ > Sj{P) \\a\\, we find 



l!«ll < 



1 



and the result follows from Lemma [2.41 D 



2.3 Complexes 

Given a finite set P, an abstract simplicial complex is a set of subsets K C 2^ such that if 
a G K, then every subset of a is also in K. The Delaunay complexes we study are abstract 
simplicial complexes, but their simplices carry a canonical geometry induced from the inclusion 
map t : P ^-> M™. (We assume l is injective on P, and so do not distinguish between P and 
t(P).) For each abstract simplex a E K, we have an associated geometric simplex conv(t(CT)), 
and normally when we write a £ K, we are referring to this geometric object. Occasionally, 
when it is convenient to emphasise a distinction, we will write l{(j) instead of cr. 

Thus we view such a iiT as a set of simplices in M™, and we refer to it as a complex, but 
it is not generally a (geometric) simplicial complex. A geometric simplicial com,plex is a finite 
collection G of simplices in M.^ such that if cr G G, then all of the faces of a also belong to 
G, and if ct, ct G G and a H a ^ 0, then r = cr n ct is a simplex and t < a and t < a. 
Observe that the simplices in a geometric simplicial complex are necessarily non-degenerate. An 
abstract simplicial complex is defined from a geometric simplicial complex in an obvious way. 
A geom,etric realization of an abstract simplicial complex X is a geometric simplicial complex 
whose associated abstract simplicial complex may be identified with K. A geometric realization 
always exists for any complex. Details can be found in algebraic topology textbooks; the book 
by Munkres |Mun84| for example. 

The dimension of a complex K is the largest dimension of the simplices in K. We say that 
K is an ?n-complex, to mean that it is of dimension m. The complex 7^ is a pure m-complex if 
it is an ?Ti-complex, and every simplex in K is the face of an 77i-simplex. 

The carrier of an abstract complex K is the underlying topological space \K\, associated with 
a geometric realization of K . Thus if G is a geometric realization of K, then \K\ = Ucrec '^- ^'^'^ 
our complexes, the inclusion map t induces a continous map t : |i^| — >■ M™, defined by barycentric 
interpolation on each simplex. If this map is injective, we say that K is embedded. In this case 
i also defines a geometric realization of K, and we may identify the carrier of K with the image 
of L. 

A subset K' <Z K is a. subcomplex of K if it is also a complex. The star of a subcomplex 
K' (Z K i& the subcomplex generated by the simplices incident to K' . I.e., it is all the simplices 
that share a face with a simplex of K' , plus all the faces of such simplices. This is a departure 
from a common usage of this same term in the topology literature. The star of K' is denoted 
star(ii'') when there is no risk of ambiguity, otherwise we also specify the parent complex, as in 
sta,v{K'; K). A simple example of the star of a complex is depicted in Figure [2J 

A triangulation of P C K™ is an embedded complex K with vertices P such that |A'| := 
conv(P). A complex K is & j -manifold complex if the star of every vertex is isomorphic to the 
star of a triangulation of MP . In order to exploit the local nature of the definition of a manifold 
complex, it is convenient to have a local notion of triangulation for the star of a vertex in K, 
even if the whole of K is not a triangulation of its vertices: 
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Figure 2: The star of a subcomplex K' C K is the subcomplex star(i^') C K that consists all 
the simplices that share a face with K' (this includes all of K' itself) , and all the faces of these 
siniplices. Here we show an embedded 2-coniplex, with all 2-siniplices shaded. The subcomplex 
K' consists of the two indicated triangles, and their faces (blue). The simplices of stsa-{K') are 
shown in bold (red and blue). The other simplices do not belong to stai-{K') (black). 



2. 
3. 

4. 



Definition 2.5 (Triangulation at a point) A complex if is a triangulation at p ^ M™ if: 
1. p is a vertex of K. 
star(p) is embedded. 
p lies in int |star(p)|. 
For all T E K, and a E star(p), if (intr) n cr 7^ 0, then t E star(p). 

In a general complex Condition |4] above is not a local property, however in the case of Delaunay 
complexes that intersts us here, local conditions are sufficient to verify the condition, as we will 
show in Section p?. 2.1 1 Observe also that Condition |4] also precludes intersections with degenerate 
simplices, since such a simplex would have a face that violates the conditon. 

If (T is a simplex with vertices in P, then any map C : P — > P C M™ defines a simplex C,{<t) 
whose verticies in P are the images of vertices of cr. If X is a complex on P, and X is a complex 
on P, then C induces a simplicial map K ^ K if C{a) E K for every a E K. We denote this map 
by the same symbol, C. We are interested in the case when C is an isomorphism, which means it 
establishes a bijection between K and K. We then say that K and K are isomorphic, and write 

i ~ . . 

K ^ K, or K ^ K if we wish to emphasise that the correspondence is given by C- 



\K\ 



K 



by barycentric inter- 



A simplicial map C, : K -^ K defines a continuous map C, 
polation on each simplex a E K . We observe the following consequence of Brouwer's invariance 
of domain: 

Lemma 2.6 Suppose X is a complex with vertices P C M™, and K a complex with vertices 
P C K™. Suppose also that K i&a. triangulation aX p E P, and that C : P — > P induces an injective 
simplicial map star(p) — > star(C(p)). If star(C(p)) is embedded, then 

C(star(p)) = star(C(p)), 
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and C(p) is an interior point of P. 

Proof We need to show that star(C(p)) C C(star(p)). Since star(p) is embedded, C defines a 
continuous map ( : |star(p)| — >• |star((^(p))| that is injective on each simplex. Since star((^(p)) is 
also embedded, this continuous map is injective on |star(p)|. Since K is a, triangulation at p, there 
is an open ball B centred at p such that B C int |star(p)|. Then C|b : R™ D S ^^ ((B) C K™ is 
a homeomorphism by Brouwer's invariance of domain |Dug66[ Ch. XVII]. It follows that ({p) is 
an interior point of P. 

Suppose a € star(C(p)) and C(p) is a vertex of a. Then, since a is not degenerate, there is a 
point X G C(^) 1^ int a, and from the above argument, x also lies in the interior of some simplex 
f e C(star(p)) C star{({p)). Since star(C(p)) is embedded, f n cr is a face of <t and of f, but since 
X is in the interior of both simplices, it must be that f = a. Thus a € ^(star(p)). 

If CT S star(C(p)), then there is some t S star(C(p)) such that ({p) is a vertex of r and a < t. 
Since t G C(star(p)), we also have a £ (^{sta,r{p)), by the definition of a simplicial map. D 



3 Parameterized genericity 



In this section we examine the Delaunay complex of P C M™ , taking the view that poorly-shaped 
simplices arise from almost degenerate configurations of points. We introduce the concept of a 
protected Delaunay ball, which leads to a parameterized definition of genericity. We then show 
that a lower bound on the protection of the maximal simplices yields a lower bound on their 
thickness. 

3.1 The Delaunay complex 

An empty ball is one that contains no point from P. 

Definition 3.1 (Delaunay complex) A Delaunay ball is a maximal empty ball. Specifically, 
B = -Br™ (x; r) is a Delaunay ball if any empty ball centred at x is contained in B. A simplex a 
is a Delaunay simplex if there exists some Delaunay ball B such that the vertices of a belong to 
dB n p. The Delaunay complex is the set of Delaunay simplices, and is denoted Del(P). 

The Delaunay complex has the combinatorial structure of an abstract simplicial complex, but 
Del(P) is embedded only when P satisfies appropriate genericity requirements, as discussed in 



Section 3.2 Otherwise, Del(P) contains degenerate simplices. We make here some observations 
that are not dependent on assumptions of genericity. 

The union of the Delaunay simplices is conv(P). A simplex a G Del(P) is a boundary simplex 
if all its vertices lie on 9conv(P). We observe 

Lemma 3.2 (Maximal simplices) If aff(P) = M™, then every Delaunay j-simplex, a, is a face 
of a Delaunay simplex a' with diniaff(cr') ~ m. In particular, if j < m, then cr is a face of a 
Delaunay ?Ti-simplex. If a is not a boundary simplex, and dim aff (cr) < m, then there are at least 
two Delaunay (j + l)-simplices that have cr as a face. 

Proof Suppose dim aff (cr) < ni. Let B = i3Rm(c; r) be a Delaunay ball for a. Let £ be the line 
through c and C{a). If c = C{(t), let £ be any line through c and orthogonal to aff(cr). There 
must be a point c Cz £ such that the circumscribing ball for a centred at c is not empty. If this 
were not the case, we would have aff(cr) = aff(P), and thus dim aff (P) < m. It follows then 
(from the continuity of the radius of the circumballs parameterized by £) , that there is a point 
c' e [c, c] that is the centre of a Delaunay ball for a simplex a' that has cr as a proper face. The 
first assertion follows. 
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Figure 3: Lemma 3.2 If the afHne hull of a is not full dimensional, then a Delaunay ball has 



freedom to expand, and a must be the face of a higher dimensional Delaunay simplex. 

The second assertion follows from the same argument, and the observation that if a is not on 
the boundary of conv(P), then there must be non-empty balls centred on £ at either side of c. If 
p (z P \ aff (cr) is on the boundary of an empty ball centred at one side of c, by the intersection 
properties of spheres, it cannot be on the boundary of an empty ball centred on the other side 
of c. Thus there must be at least two distinct Delaunay (fc + l)-simplices that share cr as a face. 

D 



Lemma |3.2| gives rise to the following observation, which plays an important role in Sec- 
tion |3.3| where we argue that protecting the Delaunay m-simplices yields a thickness bound on 
the simplices. 

Lemma 3.3 (Separation) If t e Del(P) is a j-simplex that is not a boundary simplex, and 
q (z P \ T, then there is a Delaunay rn-simplex cr™ which has r as a face, but does not include q. 

Proof Assume j < m, for otherwise there is nothing to prove. If cr = g * r is not Delaunay, 
the assertion follows from the first part of Lemma |3.2| Assume a is Delaunay and let ct'" be 
a Delaunay m-simplex that has cr as a face. Thus cr™ = q * a^"-^^ for some Delaunay (to — 1)- 
simplex, (t™~'^. Since r < a™~^ does not belong to the boundary of conv(P), neither does a"^~^ 



so by the second part of Lemma 3.2 there is another Delaunay TO,-simplex a^"' that has a™' ^ 



(and therefore r) as a face. Since a"^ is distinct from (t™, it does not have q as a vertex. D 

3.1.1 The Delaunay complex in other metrics 

We will also consider the Delaunay complex defined with respect to a metric d on M™ which 
differs from the Euclidean one. Specifically, if P C t/ C M™ and d : U x U —> M. is a, metric, then 
we define the Delaunay complex Delrf(P) with respect to the metric d. 

The definitions are exactly analogous to the Euclidean case: A Delaunay ball is a maximal 
empty ball B{x;r) in the metric d. The resulting Delaunay complex Delc;(P) consists of all the 
simplices which are circumscribed by a Delaunay ball with respect to the metric d. The simplices 
of Deld(P) are, possibly degenerate, geometric simplices in E™. As for Del(P), Delrf(P) has the 
combinatorial structure of an abstract simplicial complex, but unlike Del(P), Deld(P) may fail 
to be embedded even when there are no degenerate simplices. 
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3.2 Protection 

A Delaunay simplex a is 5-protected if it has a Delaunay ball B such that djRm {q, dB) > 5 for all 
g g P\a. We say that i? is a ^-protected Delaunay ball for a. If r < cr, then B is also a Delaunay 
ball for r, but it cannot be a (5-protected Delaunay ball for r. We say that u is protected to mean 
that it is (5-protected for some unspecified 5 >{). 




Figure 4: A Delaunay simplex a is ^-protected if it has a Delaunay ball Bk™ (c; r) such that 

BM™(c;r + (5)n(P\(T)=0. 



Definition 3.4 (J-generic) A finite set of points P C M™ is 5-generic if aff(P) = M™, and all 
the Delaunay m-simplices are (5-protected. The set P is simply generic if it is (5-generic for some 
unspecified S > 0. 



Observe that we have employed a strict inequality in the definition of (5-protection. In particular, 
a (5-generic point set is generic even when 5 — 0. In order for the quantity 5 to be meaningful. 
it should be considered with respect to a sampling radius e for P. We will always assume that 
S<e. 

In his seminal work, Delaunay |Del34| demonstrated that if there is no empty ball with m + 2 
points from P on its boundary, then Del(P) is realized as a simplicial complex in W\ In other 
words Del(P) is an embedded complex, and in fact it is a triangulation of P, the Delaunay 
triangulation. If P is generic according to Definition |3.4| then Delaunay's criterion will be 
met. This is obvious if there are no degenerate m-simplices, and Definition 3.4 ensures that a 



degenerate ?Ti-simplex cannot exist in Del(P), as shown by Lemma 3.5 below. 

In particular, if P is generic if and only if there are no Delaunay simplices with dimension 
higher than m. We can say more. There are no degenerate Delaunay simplices. This can be 



inferred directly from Delaunay's result [Del34| . but is also easily established from Lemma 3.2 



In Section 3.3 we will quantify this observation with a bound on the thickness of the Delaunay 
simplices. 

The (5-generic assumption means that all the Delaunay m-simplices are (5-protected, but the 
lower dimensional Delaunay do not necessarily enjoy this level of protection. The fact that 
there are no degenerate Delaunay simplices implies that all the simplices of all dimensions are 
(5-protected for some (5 > 0. 
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3.2.1 Local Delaunay triangulation 

Delaunay avoided boundary complications by assuming a periodic point set, but we are par- 
ticularly interested in the case where the point sets come from local patches of a well-sampled 
compact manifold without boundary. Periodic boundary conditions are not appropriate in this 
setting, but this is not a problem because, as we show here, Delaunay's argument applies locally. 
Delaunay's proof that the Delaunay complex of a generic periodic point set is a triangulation 
of M™ consists of two observations. First it is observed that if two Delaunay simplices intersect, 
then they intersect in a common face. This shows that Del(P) is embedded. The argument is 
not complicated by the presence of boundary points: 

Lemma 3.5 (Embedded star) Suppose aff(P) = M™ and p E P. If all the m-simplices in 
star(p; Del(P)) are protected, then star(p; Del(P)) is embedded, and it is a pure m-complex. 

Proof We first observe that the m-simplices in star(p) are not degenerate. If a"^ is degenerate, 
then by Lemma 3.2 there is a simplex r with afF(T) = E™, and a™' < r. We have t £ star(p), 
since p £ t. An affinely independent set of to -I- 1 vertices from r defines a non-degenerate 
?77,-simplex (j™ < t, and since its unique circumball is also a Delaunay ball for r, it cannot be 
protected, a contradiction. 

Now suppose that a,T G star(p) and cr n r 7^ 0. We need to show that they intersect in a 
common face. By Lemma |3.2[ we may assume that a and r are m-simplices. Let Bi and B2 be 
the Delaunay balls for a and r, and let H be the (to. — l)-fiat defined by a.S{dBi n 9i?2)- Since 
Bi and B2 are empty balls, if a ^ t, H separates the interiors of a and r, and thus they must 
intersect in H, i.e., at the common face defined by the vertices in dBi D dB2. □ 

The second observation Delaunay made is that, in the case of a periodic (infinite) point set, 
every {m — l)-simplex is the face of two TO,-simplices (Lemma 3.2 1. The implication here is that 
Del(P) cannot have a boundary, and therefore must cover M™. Here we flesh out the argument for 
our purposes: If an embedded finite complex contains m-simplices then its topological boundary 
must contain [m— l)-simplices. We first observe that the topological boundary of an embedded 
complex is defined by a subcomplex: 



Lemma 3.6 (Boundary complex) If K is an embedded complex in IE 

boundary of \K\ C M™ is defined by a subcomplex: d\K\ — \\)(\.{K)\, 
hA{K) <Z K is called the boundary complex of K. 



7", then the topological 
where the subcomplex 




y 



a' 



Figure 5: Diagram for the proof of Lemma 3.6 



Proof Suppose x £ d\K\. Then x E inter-' for some a^ £ K. We wish to show that a^ C d\K\. 
Suppose to the contrary that y £ inter-', but y does not belong to d\K\. This means that 

y £ iiA\K\. 
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Consider the segment £ = [x,y] C int cr^. 
siniplices that do not contain aK Let 



Let Z d K he the subcomplex consisting of those 



ri = min(iR™(£, cr). 



Since y £ mt\K\, we may 



Choosing r < ri, and x' G i?Rm(a;;r) \ \K\, let y' = y + {x' — 
assume that r is small enough so that y' G int |A'|. 

Consider the segment £' = [x',y']. By construction. £' \Z\ = 0. However, consider the point 
u e int £' that is the point in / n |iir| that is closest to x' . The point u lies in the interior of some 
simplex a' G K, but we cannot have a^ < a' . Indeed if this were the case, x' would lie in aff (cr'), 
and so u G da', contradicting the assumption that u € inta'. 

But this means that a' E Z, which contradicts the fact that £' n \Z\ = 0. Therefore we must 
have y E d\K\ for all y G int a^ . 

Finally, observe that if t < cr-', then t C d\K\, since d\K\ is closed. D 

Lemma 3.7 (Pure boundary complex) If iiT is a (finite) pure m-complex embedded in M™, 
then its boundary complex is a pure (m — l)-complex. 




Figure 6: Diagram for the proof of Lemma 3.7 



Proof Since K is finite, hd{K) is nonempty; it contains at least the vertices in dcoiw{\K\). 

We will show that if cr-' G hd{K), is a j-simplex, with < j < m—1, then there is a cr*^ g hd{K) 

with a^ < CT*^. The result then follows, since hd{K) cannot contain m-simplices, because K is 

embedded. 

Suppose a^ £ hd(K), and x E intcr^. Let Z C K he the subcomplex consisting of siniplices 

that do not contain a^ , and let 

r = min dom (x,a). 
aez 

Let B — B-Rm(x;r), and choose y E B\ \K\. Let F he the (m — j)-dimensional affine space 
orthogonal to aff(cr-') and containing y, and let 5™^-?'^! = F n 9i?Rm(a;; r'), where r' = \\x — y\\. 
See Figure [6) 
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Since K is pure, there is an m-simplex a™ with a^ < ct™. We have cr^nS'™^^^"'^ 7^ 0. Indeed, 
choose w € int cr™, and u € a^ different from x, and observe that the plane Q defined by x,w,u 
intersects B f] cr™ in a semi-disk, by construction of B. By the construction of 5*'"^^^^, it must 
intersect this semidisk. 

Let z £ S™~^~^ be a point that minimises the geodesic distance in 5''"^-'^^ to y. Then 
z e d\K\. Thus z E intCT*^ for some a'^ E bd(i^), and since z E B, a^ cannot belong to Z. Thus 
a^ < a'', but since S"""^"i n a^ = 0, we have a^ < a^. D 



From Lemma |3.2| Lemma |3.5[ and Lemma |3.7| one can verify that if P is generic then 
9conv(P) = 9|Del(P)|, and thus obtain the standard result that Del(P) is a triangulation of P. 
However, we are interested localizing the result, without the assumption that the entire point 
set is generic. We have the following local version of Delaunay's triangulation result: 

Lemma 3.8 (Local Delaunay triangulation) If p G P is an interior point, and the Delaunay 
777,-simplices incident to p are protected, then Del(P) is a triangulation at p. 



Proof By Lemma 3.5 star(p) is a pure m-complex, and it is embedded. It follows then from 



Lemma [3. 7[ that the boundary complex bd(star(p)) is a pure (ni — l)-complex. Thus p cannot 



belong to bd(star(p)). Indeed, it follows from Lemma 3.2 that any (m— l)-simplex a E star(p) is 
the face of at least two m-simplices in Del(P), and if p G a, then both of these m-simplices belong 
to star(p), and are embedded, with intersection a. Thus p cannot belong to an (m — l)-simplex 
in bd(star(p)), and therefore p E int |star(p)|. 

It remains to verify Condition |4] of Definition |2.5[ The argument is similar to the proof of 



Lenuna 3.5 Suppose x E (intr) n a for a E star(p). We may assume that a is an ?7i-simplex. 
Then consider the Delaunay balls Bi for a and B2 for r, and let H be the (m — l)-flat defined 
by afF(9J5i n dB2). Since a is protected, x can lie in the interior of a only if r = cr. Otherwise, 
X E H, and therefore all vertices of r lie in H, and r is a face of a. D 

3.2.2 Safe interior simplices 

We wish to consider the properties of Delaunay triangulations in regions which are comfortably 
in the interior of conv(P), and avoid the complications that arise as we approach the boundary 
of the point set. We introduce some terminology to facilitate this. 

If none of the vertices of a lie on 9conv(P), then it is an interior simplex. We wish to identify 
a subcomplex of the interior simplices of Del(P) consisting of those simplices whose neighbour 
simplices are also all interior simplices with small circumradius. An interior simplex near the 
boundary of conv(P) does not necessarily have its circumradius constrained by the sampling 
radius. However, we have the following: 

Lemma 3.9 If P is an e-sample set, and <j E Del(P) has a vertex p such that dR™ {p, 9conv(P)) > 
2e, then -R(cr) < e and a is an interior simplex. 

Proof Let B^m (c; r) be a Delaunay ball for a. We will show r < e. Suppose to the contrary. 
Let X be the point on [c,p] such that (iRm(j5, x) = e. Then p is the closest point in P to x, 
and so the sampling criteria imply that ds,"^{x,dconv{P)) < e. But then (iRm(p, 9conv(P)) < 
dRm (p, x) + dRm {x, 9conv(P)) < 2e, contradicting the hypothesis on p. 

Thus r < e, and it follows that a is an interior simplex because if g G cr, then ds^m (p, q) < 
2r < dM™(p,9conv(P)). D 

This suggests the following: 
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Definition 3.10 (Deep interior points) Suppose P C M™ is an e-sample set. The subset 
P/ C P consisting of all p G P with dRm(p, 9conv(P)) > 4e is the set of deep interior points. 



By Lemma [3. 9 1 all the simplices that include a deep interior point, as well as all the neighbours 
of such simplices, will have a small circumradius. For technical reasons it is inconvenient to 
demand that all the Delaunay m-simplices be (5-protected. We focus instead on a subset defined 
with respect to a set of deep interior points: 




Figure 7: If P is (5-generic for Pj, then the safe interior simplices are the simplices in star(Pj). 
Here Pj consists of the two large vertices (blue). They must be at least 4e from (9conv(P) 
(which is not depicted in the figure). The safe interior simplices are shaded. All the simplices in 
star(star(Pj)) are (5-protected. These simplices have bold outlines (red), but are not necessarily 
shaded. 



Definition 3.11 (J-generic for Pj) The set P C M™ is S-generic for Pj if Pj C Pj and all 

the m-simplices in star(star(Pj;Del(P))) are (5-protected. The safe interior simplices are the 
simplices in star(P,7; Del(P)). 

Thus the safe interior simplices are determined by our choice of Pj C P/, and our protection 
requirements ensure that all the TTi-simplices that share a face with a safe interior simplex are 
(5-protected and have a small circumradius. A schematic example is depicted in Figure [7] 



3.3 Thickness from protection 

Our goal here is to demonstrate that the safe interior simplices on a (5-generic point set are To- 
thick. If (5 = i/ge, for some constant I'o < 1, then we obtain a constant Tq which depends only 
on i/Q. The key observation is that together with Lemma |3.3[ protection imposes constraints on 
all the Delaunay simplices; they cannot be too close to being degenerate. In the particular case 
that j — 0, Lemma |3.3| immediately implies that the vertices of the safe interior simplices are 
(5-separated: 

Lemma 3.12 (Separation from protection) If P is (5-generic for Pj, then L{<t) > 5 for any 

safe interior simplex a. 
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(a) H separates c and q 



(b) q and c on same side of H 



Figure 8: Diagram for Lemma 3.13 (a) When H separates q and c then d]g^d{q,q*) > S. (b) 
Otherwise, a lower bound on the distance between q and its projection q* on H is obtained by 
an upper bound on the angle Zqab. 



Lemma 3.13 Suppose that B = Br™ (c; r) is a Delaunay ball for cr = q * r with r < e and that 
L{t) > A for some A < e. Suppose also that t < a' and that a is not a face of a'. 
If B' is a (5-protected Delaunay ball for a', and H = a.S{dB D dB'), then 

d^^{q,H)>^{X + S). 

It follows that, if P is (5-generic for Pj, with sampling radius e, and r is a safe interior simplex, 
then 

D{q,(J) > -^^■ 

Proof Let B' — BRm(c';r') be the (5-protected Delaunay ball for a'. Our geometry will be 
performed in the plane, Q, defined by c, c', and q. This plane is orthogonal to the (m— l)-flat H, 
and it follows that the distance dRm (g, H) is realized by a segment in the plane Q: the projection, 
q* , of g onto H lies in Q, and df^m(q^H) = dRm(g, g*). 

If H separates q from c, then cJi?' separates g from g*, and (iR.n(g, g*) > dg^m^q^dB') > S, 
since S' is (5-protected (Figure ^[a)[ ). The lemma then follows since A and 6 are each no larger 
than e. Thus assume that g and c lie on the same side of H, as shown in Figure ^[b) Let 
S" = Q n dB', and S ^ QndB, and let a and b be the points of intersection S'nS. Thus iJ n Q 
is the line through a and b. 

We will bound dRm(g, g*) by finding an upper bound on the angle 7 = Zqab. This is the 
same as the standard calculation for upper-bounding the angles in a triangle with bounded 
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circumradius to shortest edge ratio. Without loss of generahty, we may assume that 7 > Aqba, 
and we will assume that 7 > 7r/2 since otherwise g* £ B' and thus djf™ (<Z; 9*) > ^ S'lid the lemma 
is again trivially satisified. 

Since d-srn{a,q) > S, we have dRm{q,q*) = dRm(a, g)sin7 > 5 sin 7. Also note that d{a,b) > 
2R{t) > L{t) > A. Let a = Zqac. Then cos a = "'"2^°''^ > ^, which means that a < 
arccos ^ < f • Similarly, with /? — Zcab, we have /3 < arccos ^j < f • Thus | <7 = a + /3<7', 
where 

5 A 

7 = arccos h arccos — . 

Since sin 7 > sin 7', when f < 7 < 7' < tt, we have 

dR^{q,q*) > <5sin7 > (5 sin 7' 

^ ■ f S ^ 

= sm arccos h arccos — 

V 2e 2e 

^.^^•^ '^\ S . f A 

> — sm arccos — H sm arccos — 

- \2e V 2ey 2e V 2e 

where the last inequality follows from A < e and (5 < e. 

Since afF(T) C -ff , it follows that I?(g, u) > dij™ (g, iJ), and if P is 5-generic for Pj, then A > (5, 
and Lemma |3.3| ensures that there is a (5-protected a' that contains r but not q. D 



We thus obtai n a bo und on the thickness of the safe interior simplices when P is (5-generic for 
P 7 . Since Lemma 



3.13 



yields a lower bound of ^^^ on the altitudes of the safe interior simplices, 



and since A(cr) < 2e, we have that T(cr) > ^^2 for all safe interior a. If (5 = v^e, we obtain a 
constant thickness bound. 

Theorem 3.14 (Thickness from protection) If P C M™ is (5-generic for Pj with 6 = t'oe, 
where e is a sampling radius for P, then the safe interior simplices are Tp-thick, with 

JLo — — -, — • 



4 Delaunay stability 

We find upper bounds on the magnitude of a perturbation for which a protected Delaunay ball 
remains a Delaunay ball. We consider both perturbations of the sample points in Euclidean 
space, and perturbations of the metric itself. The primary technical challenge is bounding the 
effect of a perturbation on the circumcentre of an TTi-simplex. We then find the relationship 
between the perturbation parameter p and the protection parameter (5 which ensures that a 
(5-protected Delaunay simplex will remain a Delaunay simplex. 

4.1 Perturbations and circumcentres 

As expected, a bound on the displacement of the circumcentre requires a bound on the thickness 
of the simplex. 
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4.1.1 Almost circumcentres 

If a is thick, a point whose distances to the vertices of a are all almost the same, will lie close to 
Nia). 

Lemma 4.1 If ct = [po, ■■■,Pk] C M'" is a non-degenerate /c-simplex, and x G M™ is such that 

<e foralH,je[0,...,A;], (2) 



\Pi -xW'^ - Wpj -xf 



then there is a point c E N{a) such that \\c — x\\ < rj, where 

e 



V 



2T((t)A((t)' 



In particular, if a is an m-simplex then x S B]^m [C{(t); rf). 

If the inequalities in Equations Q are made strict, then the conclusions may also be stated 
with strict inequalities. 

Proof First suppose k ^ m. The circumcentre of a is given by the linear equations || C(cr) — Pi\\ 



ip,-p,fCia)^l{\\p£-\\pof). 



|lC(a)-poll ,or 

( 

Letting b be the vector whose i^^ component is defined by the right hand side of the equation, 
and letting P be the m x m matrix, whose i*'^ column is {pi — po), we write the equations in 
matrix form: 

P'^C{a) = b. (3) 

Without loss of generality, assume pq is the vertex that minimizes the distance to x. Then, 



defining Xa to be the vector whose i^^ component is |(||pi 
IIpo ^ 2;|| + 2(xa)i, and we find 



l|Po-a;|| ), we have \\p^ - x\\ 



P X = b — x„ 



(4) 



From Equations (l3| and Q we have 



\\C{a)-x\\^ (PT) '^^ 



< P' 



From Equation (pi, it follows that \\xa\ 



< 



/^r 



and from Lemmas 



2.3 



and 



2.4 



we have 1 1 P "'" 1 1 < 



(•y/mT(cr)A((T)) , and thus the result holds for full dimensional simplices. 

If (J is a /c-simplex with k < m, then we consider x, the orthogonal projection of x into aS{a). 
We observe that x also must satisfy Equation (|2l, and we conclude from the above argument 
that ||C(cr) — £'11 < rj. Then letting c = C{a) + {x — x) we have that c £ N{a) and \\c — x\\ < r). 

D 

It will be convenient to have a name for a point that is almost equidistant to the vertices of 
a simplex: 

Definition 4.2 A (^-centre for a simplex a = [po, . . . ,pk] C M™ is a point x that satisfies 

\\pi — x\\ — \\pj ~ x\\\ < ^ for all i,j < k. (5) 



With a bound on the distance from x to the vertices of a, Lemma |4.1| can be transformed into a 
bound on the distance from a ^-centre to the closest true centre in N{a): 
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Lemma 4.3 li a — [po^ ■ ■ ■ ,Pfe] C M™ is non-degenerate, and for some <^ > the point x € 
is a ^-centre such that 

\\pi — a;|| < e for all i,j < fc, 

then there exists a c G N{(j) such that ||a; — c\\ < rj, where 

'' T(a)A(a)- 
In particular, if a is an 771-siniplex, then x G B^m(^C{(7);ri). 
Proof Let R — maxi ||pi — x\\ and r = miui \\pi — a;||. Then 

i?2 - r^ = (i? + r)(i? - r) < 2?(i? - r) < Tti, 
and the result then follows from Lemma |4. II 



n 



4.1.2 Circumcentres and metric perturbations 

We will show here that for an Tg-thick m-simplex a in K™, and a metric d that is close to dR™, 
there will be a metric circumcentre c near C((7). We require the metric d to be continuous in 
the topology defined by d^-m. Henceforth, whenever we refer to a perturbation of the Euclidean 
metric, this topological compatibility will always be assumed. 




Figure 9: The maps /e and / (described in the main text) map circumcentres to the origin. Since 
/~^(0) contains a unique point, and fe{dB) lies far from the origin, a consideration of the degree 
of the mappings, together with the fact that fe{dB) and f{dB) are close, reveals that /^^(O) 
cannot be empty, and thus B must contain a circumcentre of tr'" with respect to the metric d. 

The proof is a topological argument based on considering a mapping into M™ of a small ball 
around the circumcentre of a. The mapping is based on the metric and is such that circumcentres 
get mapped to the origin. In the mapping associated to the Euclidean metric, points that get 
mapped close to the origin are ^-centres, and since the ^-centres are in the interior of the ball, 
the boundary of the ball does not get mapped close to the origin. A small perturbation of 
the metric yields a small perturbation in the mapping, and so we can argue that there is a 
homotopy between the mapping associated with the Euclidean metric and the one associated 
to the perturbed metric, such that no point on the boundary of the ball ever gets mapped to 
the origin. The situation is depicted schematically in Figure [9] A consideration of the degree of 
the mapping allows us to conclude that the ball must contain a circumcentre for the perturbed 
metric. 

We will demonstrate the following: 
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Lemma 4.4 (Circumcentres: metric perturbation) Let U C M™, and let d : U x U ^ M. 
be a continuous metric with respect to the topology defined by djjm , and such that for any 
x,y £U with (iRm(x,j/) < 2e, we have \d{x,y) — d^m{x,y)\ < p, with 

If cr = [poi ■ • ■ iPm] C [/ is an Tg-thick m-simplex with R{<t) < e, and L{a) > fiae, and such that 
c'r™ (Pjj f^t^) ^ 2e. then there is a point 



c 



e S = i3R"^ (C(cr) ; 77) with i] = 



ToMo' 

and such that d{c,Pi) ~ d{c,pj) for all Pi^Pj G a. 

In order to prove Lenuna |4.4[ we will use a particular case of Lemma |4.3[ 

Lemma 4.5 Suppose <j is an Tg-thick ?7i-simplex such that L{a) > fiQe. If a; is a ^-centre for a 
with d-Rm(^x,p) < 2e for all p ^ cr, then x G i?M™ (C(cr); rj), where 77 = -^ ^ . 

Let i? = i?Km(C(cr);7y) be the open ball which contains the ^-centres for a. We will show that if 
^ = 4p, then a circumcentre c for cr with respect to d will also lie in B. However, we make no 
claim that c is unique. Note that B C U. 

Consider the function fe'-B^ M™ given by 

Je(.x) ^ {dK,^[x,pi)~ d^^{x,po),...,d^rr.{x,pm)- dwL^[x,po)) . (6) 

Observe that /g maps the circumcentre of a, and only the circumcentre, to the origin: /^^(O) = 
{C{a)}. 

We construct a similar function for the metric d, 

f{x) ^ {d{x,pi)-d{x,po),...,d{x,p„i)-d{x,po)) , (7) 

and we will show that there must be a c G /^^(O) C B. We first show that there is a homotopy 
between / and /e such that the image of dB never touches the origin: 



Lemma 4.6 Under the hypotheses of Lemma 4.4 if C = 4/0 < "^°^ , then there is a homotopy 



Tafige 



F : B X [0, 1] -^ E™ between f^{x) = F{x, 0) and f{x) = F{x, 1) with F{x, t) ^ for all x e dB 
and te [0,1]. 

Proof We define the homotopy F :Bx[0,l]^ M'" by F(a;, i) = (1 - t)fe{x) + tf{x). By the 
bounds on ^ and R{(j), for every x £ B, and p £ a, we have dR»».(a;,p) < y"^ — hi?(o') < 2e. Thus 
it follows from Lemma 4.5 that x € 95 cannot be a ^-centre. 

It is convenient to consider the max norm on M™ defined by the largest magnitude of the 
components: ||/e(a;)||oQ = max^ \fe{x)i\. (This gives us a better bound than working with the 
standard Euclidean norm.) If ||/e(y)|loo < f' then y must be a ^-centre. Indeed, we would have 
llbi-y|| - lbi-y|!l < llb»-y|| - lbo-2/||| + ll|Po-2/|| - \\Pj-y\\\ < f + f =Cfor alH and j. 
Thus, since x £ dB is not a ^-centre, we must have ||/e(a;)||g^ > |. 

Also, from the hypothesis on d, we have \\feix) — fix)\\^ < 2/? = |, for all x G dB. It follows 
that \\Fix,t)\\^ > ||/e(x)||^ - t \\f{x) - fe{x)\\^ > when x e OB. O 

We will need the following observation: 
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Lemma 4.7 The origin is a regular value for the function /e defined in Equation (Ig]). 

Proof Choose a coordinate system such that C{a) E B is the origin. We show by a direct 
calculation that det J(/e)o y^ 0, where J(/e)o is the Jacobian matrix representing the derivative 
of /e in our coordinate system. 

Let Pi = {pio, ...,pim) for all pi e {po, ...,pm}. For x = {xi, . . . , a;„) G M", let 
fe{x) = {fo{x), ■■■, fm{x)) , where 



We find 



and thus 



Mx) = WPz -x\\ - Wpo ~x\\ 



dxj 



Afh 



\ k=l 



POj - PiJ 

R{a) ' 



R{ay 



\ k=l 



(8) 



where as usual P is the matrix whose columns are pi —po- Since vol(o'™) — ' ° \ ' , Equation 
implies 

m! vo^cr™) 



|detJ(/)o| 



R{ay 



Thus since /^ ^(0) — {0}, is a regular value for /g provided a is non-degenerate. 



n 



Lemma |4.4| now follows from a consideration of the degree of the mappings / and /g relative 
to zero. The degree of a smooth map g : B ^ M™ at a regular point p £ g{B) is defined by 

deg(g,p,S)= ^ signdet J(5);:,, 



where J{g)x is the Jacobian matrix of g at x. The exposition by Dancer [DanOO] is a good 
reference for the degree of maps from manifolds with boundary. It is shown that the definition of 
deg((7,p, i?) extends to continuous maps g that are not necessarily differentiable. li h : B -^ M™ 
is homotopic to g by a homotopy 7f : i? x [0, 1] ^ M'" such that H{x, t) ^ p for all t E [0, 1], and 
X E dB, then deg{g,p,B) = deg{h,p,B). 

Since /~^(0) — {C(cr)}, it follows from Lemma 4.7 that deg(/e, 0, B) = ±1. Then Lemma 4.6 
implies deg(/, 0,i?) = deg(/e, 0, i?), and since this is nonzero, it must be that /^^(O) 7^ 0. The 
demonstration of Lemma |4.4| is complete. 



4.1.3 Circumcentres and point perturbations 

The exact same argument as was used to demonstrate Lemma |4.4| can be used to show that 
an TO-simplex a — [po, . . . ,Pm] whose vertices lie close to a thick m-simplex a, will also have a 
circumcentre that lies close to C{a). We replace the function / defined in Equation Q by the 
function 

f{x) = (dR™(a;,pi) - dT!.^{x,po), . ■ . , dMm{x,Prn) -dRm{x,po)) , 

and the same argument goes through. We obtain: 
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Lemma 4.8 (Circumcentres: point perturbation) Suppose that a = [po, ■ ■ ■ ,Pm] is an To- 
thick ?7i-siniplcx with i?(cr) < e and L{a) > fiQe. Suppose also that a = [po, . . . ,Pm] is such that 
||Pi ~ Pill — P fo'^ all i G [0, . . . , m]. If 

P<^^, then dRn.{C{a),Cia)) < ^^ 



ToMo 

4.2 Perturbations and protection 

Suppose C : P — > P is a /7-perturbation. If cr is a (5-protected ?7i-siniplex in Del(P), then we 
want an upper bound on p that will ensure that a — Clc) is protected in Del(P). The following 
definition will be convenient: 

Definition 4.9 (Secure simplex) A simplex a <E Del(P) is secure if it is a (5-protected m- 
siniplex that is To-thick and satisfies R{cr) < e and i(cr) > /io£- 

Our stability results apply to subconiplexes of secure siniplices. the definition of which employs 
multiple parameters. For safe interior siniplices Lemma |3.12| and Theorem |3.14| allow us to 
consolidate some of these parameters with the ratio S/e: 

Lemma 4.10 (Safe interior simplices are secure) If P satisfies a sampling radius e and is 
J-generic for Pj, with 6 = vqe, then the safe interior m-simplices are secure, with fiQ = i^q, and 

iO — 4 • 

Lemma 4.11 (Protection and point perturbation) Suppose that P C M™ and <t E Del(P) 
is secure. If (^ : P — >■ P is a p-perturbation with 

, Tq^o . 

then C(cr) ~ a E Del(P) and has a ((5 — x^^ p)-protected Delaunay ball. 

Proof Let B = i3Rm(c;r) be the ^-protected Delaunay ball for a E Del(P). and let B = 
BRm{c;f) be the circumball for the corresponding perturbed simplex a. We wish to establish a 
bound on p that will ensure that B is protected with respect to P. 

Let q E P he & point not in a. We need to ensure that the corresponding q lies outside the 
closure of B, i.e., that d]^m{q,c) > r. 



Since S < e, the hypothesis of Lemma 4.8 is satisfied by p, and we have dum (c, c) < rjp, where 
77 = ^ . Thus for p E cr and corresponding p E a we have 

f < dm^ {c,p) + dR". (c, c) -I- dRm {p, p) 
< r + (77 + l)p. 

Also 

^R'" (9, c) > c^r™ (q, c) - dR-r. (5, c) - cIr™ (g, q) 
> r + 5- p{-q + l). 

Therefore q will be outside of the closure of B provided r + (5 — p[r] + 1) > r + (1 + ri)p, i.e., when 
5 > 2(77 + l)p. The result follows from the definition of 77 and the observation that po and To 
are each no larger than one. D 
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A similar argument yields a bound on the metric perturbation that will ensure the Delaunay 
balls for the m-simplices remain protected: 

Lemma 4.12 (Protection and metric perturbation) Suppose U C M'" contains conv(P) 
and (i:[/xJ7— ^Misa metric such that jdK™ {x^ y) — d{x, y)\ < p for all x,y ^ U. Suppose also 
that (J e Del(P) is secure. If 

p < ^ 



20 

ToMo' 



and dR™(p, dU) > 2e for every vertexp e a, then a also belongs to De\d{P), and has a (J— -f^^^ /o)- 
protected Delaunay ball in the metric d. 



Proof Let B = B]^m{c; r) be the Euclidean (5-protected Delaunay ball for a € Del(P), and let 
B = i?Rm (c; f) be a circumball for a in the metric d. We wish to establish a bound on p that 
will ensure that B is protected with respect to d. 

Let q £ P be a point not in a. We need to ensure that d{q, c) > f. Since S < e, the hypothesis 

8 
To Mo ' 



ensures that p < — f— , and so Lemma 4.4 yields (iRm(c, c) < jjp, where i] — y^ — . Thus for p G a 



f < d{c,p) + d{c,c) 
<{r + p) + {rip + p) 
= r + {ri + 2)p, 



and 



d{q, c) > d{q, c) — d{c, c) 
>r + 5- {■q + 2)p. 

Thus q will be outside of the closure of B provided r + 5 — {rj + 2)p > r + {rj + 2)p, i.e., when 

S>2{r] + 2)p. 

The result follows from the definition of rj and the observation that pg and To are each no larger 
than one. D 

4.3 Perturbations and Delaunay stability 

The results of Section [L2| translate into stability results for Delaunay triangulations. In the case 
of point perturbations in Euclidean space, the connectivity of the Delaunay triangulation cannot 
change as long as the simplices corresponding to the initial m-simplices remain protected. This is 
a direct consequence of Delaunay 's original result |Del34) . but we explicitly lay out the argument. 
In the case of metric perturbation, we can no longer take for granted that the Delaunay 
complex cannot change its connectivity if the m-simplices remain protected. This is because 
we are no longer guaranteed that the Delaunay complex will be a triangulation. Using the 
consequences of the point-perturbation result, we establish bounds that ensure that the Delaunay 
complex in the perturbed metric will be the same as the original Delaunay triangulation. 

4.3.1 Point perturbations 

A consequence of Delaunay's triangulation result is that if a perturbation does not destroy any 
TO-simplices in the Delaunay complex of a generic point set, then no new simplices are created 
either, and the complex is unchanged. More precisely we have: 
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Lemma 4.13 Suppose P C M™ is a generic sample set, and Q C P is a subset of interior points. 
If C : P ^- P is a perturbation such that ^(star(Q; Del(P))) C star(C(Q); Del(P)), and every 
m-simplex ct'" G C(star(Q)) is protected in Del(P). then C(star(Q)) = star(C(Q)). 



Proof Let p £ Q By Lemma 3.5 star(C(p)) is embedded, and by Lemma 3.8 Del(P) is a 
triangulation at p. Since C : P — >■ P is injective, it follows that the simplical map induced by C 
must be injective, and the result follows from Lemma |2.6| D 



Lemma |4. 1 1 1 establishes bounds on a /j-perturbation C : P ^ P which will guarant ee th at if 

also 



4.11 



Q C P, and the simplices in star(Q) are secure, then C(star(Q)) C Dcl(P). Lemma 
guarantees that, if p is small enough, the m-simplices in star(C(Q); Del(P)) will be protected 
Thus if Q consists only of interior points of P, Lemma 4.13| applies. We have the following 
stability theorem for protected Delaunay triangulations: 

Theorem 4.14 (Stability under point perturbation) Suppose P C M'" and Q C P is a 

subset of interior points such that every 77i-simplex in star(Q) is secure. If ^ : P — > P is a 
p-perturbation, with 

p < d 

^ - 18 

then 

star(Q; Del(P)) = star(C(Q); Del(P)). 

The p-relaxed Delaunay complex for P was defined by de Silva tdS08' by the criterion that 
a S Del''(P) if and only if there is a ball B = Bs.^ (c; r) such that a C B, and dftm (c, q) > r ~ p 
for all (7 e P. Thus the simplices in DeF(P) all have "almost empty", balls centred on a p-centre 
for a. We have the following consequence of Theorem |4.14| 

Corollary 4.15 (Stability under relaxation) Suppose P C M™ and Q C P is a set of interior 
points such that every 77i-simplex in star(Q) is secure. If 

, Topo . 

then 

star(Q; DeF(P)) = star(Q; Del(P)). 

Proof Suppose that a S star(Q; DeF(P)). Then there is a ball B enclosing a such that any 
point q £ B is within a distance p from dB. Project all such points radially out to OB. Th en we 
have a p-perturbation ^ : P — >■ P, and a has become a G star(C(Q);Dcl(P)). By Theorem 4.14 

star(C(Q); Del(P)) ^ star(Q; Del(P)), and therefore a € star(Q; Del(P)). D 

4.3.2 Metric perturbation 

For a perturbation of the metric, we can exploit the stability results obtained for perturbations 
of points in the Euclidean metric to ensure that no simplices can appear in star(Q; Delci(P)) that 
do not already exist in star(Q;Del(P)). 

Lemma 4.16 Suppose conv(P) C t/ c M™ and d : UxU — > Mis such that \d{x,y) — d]g,m{x,y)\ < 
p for all x,y S U. Suppose also that Q C P is a set of interior points such that every m-simplex 
a e star(Q) is secure and satisfies d^m (p, dU) > 2e for every vertex p £ a. If 

, Topo . 
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then 

star(Q;Del<i(P)) C star(Q;Del(P)). 

Proof Let B{c;r) be a Delaunay ball for simplex a S star(Q;Deld(P)). Then d{c,p) < d{c,q) 
for all p ^ (J, and g e P. By the hypothesis on d, this implies that d^m (c, p) < dffi™ (c, q) + 2p for 

D 



all p € cr and q G P, and therefore cr € Del^''(P). The result now follows from Corollary 



4.15 



The perturbation bounds required by Lemma [4.16| also satisfy the requirements of Lemma [4.12| 
This gives us the reverse inclusion, and thus we can quantify the stability under metric pertur- 
bation for subcomplexes of secure simplices in Delaunay triangulations: 

Theorem 4.17 (Stability under metric perturbation) Suppose conv(P) C [/ C M™ and 
d : C/ X [/ ^ M is such that \d{x , y) — d^m [x , y)\ < p for all x,y £ U. Suppose also that 
Q C P is a set of interior points such that every rn-simplex a £ star(Q) is secure and satisfies 
dRm (p, dU) > 2e for every vertex p & a. If 

'^- 36 ' 

then 

star(Q;Del<i(P)) = star(Q;Dcl(P)). 

Using Lemma [4.10[ and recognizing that the safe interior simplices also satisfy the distance from 
boundary requirement of Theorem |4.17| we can restate this metric perturbation stability result 
for Delaunay triangulations on (5-generic point sets: 

Corollary 4.18 (Stability under metric perturbation) Suppose P is (5-generic for Pj. with 
sampling radius e and 6 = voe. Suppose also that conv(P) C U, and rf : [/ x C/ — > M is such that 
\d{x,y) - dm^{x,y)\ < p for all x,y eU. If 



,3 ,,4 



<'^s "• 



°e, 



84 84 

then 

star(Pj;Del<i(P)) = star(Pj;Del(P)). 

5 Conclusions 

We have quantified the close relationship between the genericity of a point set, the quality of the 
simplices in the Delaunay complex, and the stability of the Delaunay complex under perturbation. 
The problem of poorly shaped simplices in a higher dimensional Delaunay complex can be seen 
as a manifestation of point sets that are close to being degenerate. The introduction of thickness 
as a geometric quality measure for simplices facilitated the stability calculations, which develop 
around a consideration of the circumcentres of a simplex in the presence of a perturbation. 

We considered a point set P C K™ meeting a sampling radius e and showed a constant bound 
on the thickness of the Delaunay simplices provided P is (5-generic with S = vqc for some constant 
uq. The question then arises: What is the least upper bound on the feasible vq as a function of 
the dimension ml 

In a companion paper |BDG13| . we develop a perturbation algorithm which produces a S- 
generic point set from a given e-sample set. In future work we will extend this to an algorithm for 
triangulating manifolds that will exploit only the local intrinsic metric properties of the manifold, 
with no requirement that it be embedded in an ambient space. 
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An important aspect of the current work that facihtates this goal, is that the triangulation 
results are localised. Since a manifold can be locally well approximated by Euclidean space, the 
objective is to fit together local Euclidean Delaunay patches where the Euclidean metric varies 
slightly between patches. This is where the stability of the Delaunay patches is important. In this 
setting we can also accommodate small variations in the sampling radius between neighbouring 
patches. Thus the algorithm will be able to triangulate sample sets whose sampling radius is 
defined by a Lipschitz sizing function. 
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